Interface resistance in ferromagnet/superconductor junctions 
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, Results of theoretical study of spin-polarized tunneling in ferromagnet/superconductor junctions 

' are presented. Spin and charge currents are calculated as a function of applied voltage and spin 

Q\ \ polarization in a ferromagnet. The model takes into account the splitting of different spin subbands 

t-H . in a ferromagnet and impurity scattering in the contact. The excess resistance of an FS contact 

i-H ' due to the charge-imbalance in a superconductor is calculated for the first time. The results have 

implications for spin-coupled magnetoresistance in ferromagnet/superconductor contacts and for 
measuring spin polarization in ferromagnets. 
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I. INTRODUCTION 

c : 

The rapidly emerging field of spin polarized transport is based on the ability of a ferromagnetic metal to conduct 
i and accumulate spin-polarized currentsoLl. Spin polarized transport between ferromagnets (F) and superconductors 
(S) received considerable attention recently because of new physical phenomena and potential device applications. 
An introduction of the hybrid structures based on a combination of ferromagnetic and superconducting materials are 
not only interesting from a fundamental point of view but can bring further advantages for devicesou. In particular, 
, t \ spin accumulation effects in superconductors may play an important role because of a number of reasons. First, 
■ due to the gap in the. excitation spectrum, the spin diffusion length in a superconductor can become quite long 
at low temperaturesacl. Second, spin accumulation can take place at a FS interface rjjJjice spin-polarized current in 
a ferromagnet has to be transformed into spinless supercurrent in a superconductomu. An important step in the 
quantitative analysis of spin accumulation and spin injection in superconductors is the knowledge of the dependence 
of the resistance of a FS interface on the spin polarization in a ferromagnet. 
£j ' Furthermore, it was recently demonstrated that a combination of F and S metals can be usee -. advantageously 
for measuring spin polarization in metallic fcrjx>magnets either by measuring T c of FS multilayersu, or by directly 
measuring the resistance of FS point contacia3'EI. So far, theoretical studies of the FS contact resistance were limited 
by calculations for ballistic FS contactsETQ It was argued inE£l that the effects of impurity scattering are quite 
important in spin-polarized tunneling since the degree of spin polarization is defined differently in ballistic and 
Q*^ ' diffusive contacts. However, no quantitative calculations on the effects of disorder have been published up to now. 
t-H , Moreover, the contribution of the contact resistance is most easily measured only in a point contact geometry, while 
for larger area contacts the contribution of an interface becomes rather small. The purpose of the present paper 
is twofold. First, we extend the theory in order to include the effects of impurity scattering in a contact. Second, 
we argue that an additional sensitive probe for spin polarization is the excess resistance R ex of a FS contact. This 
resistance is due to penetration of an electric field into a superconductor over macroscopically large charge-imbalance 
relaxation length Xq and may exceed the direct interface resistance. We show that the magnitude of R ex is sensitive 
' to spin polarization in a ferromagnet and provide an estimate for this effect. 
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We start from the derivation of the general expression for the conductance of a FS contact in the absence of impurity 
scattering (ballistic case). We consider the atomically sharp interface barrier at x = separating F metal (x < 0) and 
S metal (x > 0), modeled by a potential U(r) — H5(x) and arbitrary relation between Fermi velocities in F and S, v F ^, 
v FL and v Fs . Here H is the barrier strength parameter, v F ],i = y / 2E F ^ i /m = y / 2E F (l ± h)/m, v Fs 
where E ex , E F and E Fs are respectively exchange energy in a ferromagnet and Fermi energies in F and S metals, the 
indices |, J. refer to the spin subbands and h = E ex /E F denotes the dimensionless spin polarization in a ferromagnet. 
We assume that the effective electronic masses m F , ms are equal to the free electron mass m e , the mean free path 
is larger than the size of the contact and the pair potential is approximated by the step function A(x) = A(T)9(x), 
A(T) being the bulk pair potential in a superconductor. !_. 

Charge and spin currents can be calculated within the framework of the BTK approachEj, i.e. considering explic- 
itly Andreev and normal reflections at the FS interface and taking into account that an incoming electron and an 



Andreev reflected hole occupy opposite spin subbandalO. The electron- and hole-like excitations are represented by 
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two-component wave functions, which obey the Bogoliubov de-Gennes equations. An electron, incoming from the 
ferromagnet F into the superconductor S, is described by a plane wave with a wave vector 

i>vnc = ^ ) e ' k ^ X - W 

Due to the four-fold degeneracy of an excitation in a superconducting state, the electron is partially reflected into F 
as an electron with the opposite wave vector —k^ or as a hole 

Aefl = « ( J ) Z KlX + b ( J ) ( 2 ) 

and partially transmitted into S without branch crossing kf x , or with branch crossing k~ x 

^ rans =c^iy^*+d^ v u y-^*. (3 ) 

Here k% , fc^ and kf- are the projections of the Fermi wave vectors in two spin sibbands and in a superconductor on 
the direction x normal to the contact plane, index + (— ) refers to electron- or hole-like quasiparticles. 

The amplitudes a,b,c and d have to be determined from the matching conditions for *5fp = ipi nc + iprefi and 
^ s — Straus at the interface, x — 

vMO) = *s(0), (4) 
d^ F {0) d* s (0) 



dx dx 



= 2m e H^ F (0)/h. 



Using these conditions we find the Andreev and normal reflection coefficients, A = |a| 2 and B = |&| 2 , and the trans- 
mission coefficients with or without branch crossing, C = |c| 2 and D — |<i| 2 , respectively, which determine charge and 
spin currents in a ballistic FS contact. 

Charge current in a FS contact is given by 



1 2nh^ P] - 1 J (2ir) 

1,1 



d 2 fc,| r +oc 



l + ^i u (B)-B u (£) 



[f(E + eV)-f(E)}dE, (5) 



where Pf± — (Ep ± E ex )/2Ep = (1 ± h)/2, wf x t i i s the projection of the Fermi velocity on the direction x, f(E) 
is the Fermi distribution function, and fey is the component of the Fermi momentum parallel to the junction plane, 
which is conserved for each individual scattering process. The ratio kp\i /kp-tA provides the normalization of the 
total probability current, taking into account that Andreev scattering involves different spin subbands. 

In the limit of low temperatures T <C T c , we arrive the following expression for the charge conductance of a ballistic 
FS contact at the subgap bias voltage eV < A(T) 

r or r^r ( 1 + a2 ) p i( v n + v fi)/vfs 

(i-nn) + a (i + nn) 

Here Go = e 2 k 2 F S / Air 2 Ti is the normal state (Sharvin) conductance of the contact, S is the contact area, Tf and Tj 
are the transmission probabilities for scattering from the spin up(down) subband into a superconductor 

= 4t; FT t; Fs = 4:V Fl vp s 



r u = y/l-T hl , Z = H/hvF S , a = y/A 2 (T) - (eVf/eV. (8) 

In relevant limits the expressions (||)-(|[) agree with the results derived irJi~0Eil0, while the advantage of the 
representation (||)-((8|) is, that the charge conductance is directly expressed in terms of the individual probabilities ^ 
and therefore is particularly suitable for consideration of the impurity scattering, as explained in the next section. 
It follows from eq.(|) that in a NS contact (T T = T l = T ) at zero bias the BTK result G NS = 2G T 2 /(2 - T) 2 is 
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recovered, which yields the conductance doubling Gns — 2Go for a fully transmissive contact, T = 1. This conductance 
enhancement is suppressed in a FS contact when TV ^ Tj due to spin polarization in a ferromagnet. 

It is straightforward to extend the expressions ©-(§) to the regime of high bias voltage, eV > A(T). The results 
of numerical calculations of the dependence of charge current and zero-bias conductance are presented in Figs. 1-3 
for various values of spin polarization. It is seen that the zero-bias charge conductance is quite sensitive to the spin 
polarization h = E ex /Ep. Fig. 3 shows that for small values of h this dependence is linear, which reflects the simple 
fact that the number of transmitted electronic modes scales like 1 — h due to spin reversal by the Andreev reflection. At 
h close to 1 a nonlinearity appears in the Gfs y s h dependence due to the dependence of the transmission coefficients 
/•.. on /, fljj). 

Spin current in a FS contact is given by 



e 
2~7rft 



r,4 



g -Fl,T 



A hl (E)-B u (E) 



[f(E + eV)-f(E)] dE 



(9) 



Difference in the sign of the contributions of the Andreev coefficients A^(E) in eqs. (||) and ([)]) reflects the fact that, 
while an Andreev reflected hole carries charge in the same direction as an incoming electron, it carries spin in the 
opposite direction. 

The low temperature spin current vanishes at subgap voltages, while at eV > A(T) the spin conductance Gp S is 
given by the following expression 



u fs — 



4G /3 



[l + f3-r iri (l-(3)Y 
13 = ^/(eV) 2 - A 2 (T)/eV. 



TfV F f T l P- l v Fl 



2v Fs 



2v Fs 



(10) 
(11) 



From eq. ( |lC)| ) one can calculate numerically the spin conductance as a function of bias voltage and spin polarization. 
Figs.l, 2 show the results of numerical calculations of the dependences of the low temperature spin conductance on 
the polarization in a ferromagnet for two different values of the barrier strength parameter, Z — and Z = 1. 



III. DIFFUSIVE FS CONTACT 



In the previous section the case of a ballistic FS contact was considered, when the contact size is smaller than the 
electronic mean free path. However the latter condition is not always fulfilled in experiments, and it is therefore of 
interest to evaluate the effect of impurity scattering in a contact. As a model for a FS contact we consider two bulk 
reservoirs (S and F), which in addition to the interface potential U(r) = H5(x) are separated by the scattering region 
(a diffusive conductor) with a size smaller than the electronic mean free path. 

The expressions derived above are particularly suitable for the introduction of impurity scattering, since they allow 
straightforward application of the scattering theory. According to this theory (seec3~li3 and references therein) any 
diffusive conductor having size larger than the electronic mean free path is characterized by universal distribution of 
transmission eigenvalues t over different channels. An average conductance of a diffusive metal is then given as a sum 
of contributions of those channels, each having the conductance Go = e 2 /2nT% 

Go = ^ f 9a(t)p(t)dt, (12) 

^0 Jo 

where a =f , J, is the spin index, g a is the conductance of a channel with transmission coefficient t, G^a = e 2 N rT (Q)D a 
is the normal state conductance per spin direction, AT ff is the density of states at the Fermi level. The expression ( p^ ) 
is valid when the impurity scattering does not mix different spin directions. ■— ■ 
Function p(t) is the universal distribution function of transmission eigenvalues for different channels given byll3 

*> " 57b (13) 



and does not depend on microscopic parameters of a diffusive conductor. Eq.([13|) shows that the transmission eigen- 
values have a bimodal distribution with a peak at unit transmission and a peak at exponentially small transmission. 

As a model for the diffusive SF contact we consider two scattering regions in series: an incoming electron is first 
transmitted through the diffusive region with probability t, then crosses the FS interface with probability Tf. In 
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turn, an Andreev-reflected hole is first scattered by the interface (probability TjJ, then by the diffusive region. The 
probabilities of these two-step processes T± t 2 are given by the expression 



t + T td -tT- 



(14) 



which follows from averaging over transmission resonances between two scattering regions, assuming that all relevant 
distances exceed the electronic wave-length. 

The charge conductance in a diffusive FS contact is given by the expression eq.(^) in which the probabilities Tj.j 
should be substituted by the probabilities Ti j2 of the two-step scattering processes described above. Here we present 
the result for low temperatures and eV < A(T) 



Gfs — G 



N 



1 T 1 T 2 (l + a 2 )P l (v 2 n /v 2 Fs + l) dt 



(1 - rir 2 ) 2 +Q 2 (1 + rir 2 ) 2 tyT^T 



(15) 



where the probabilities Ti j2 are given by eq. ([l4|) , ri j2 = ^/l — Ti, 2 , a — ^/A 2 (T) — (eV) 2 /eV and Gn 
2e 2 N(Ep)D(EF) is the conductance of a contact in the unpolarized state. 

In the NS case with a transparent interface (Z=0) and vfi = vfi = ^Fs the above expression at V = yields 



G 



NS 



(V = 0) = G N f 
Jo 



tdt 



(2-t) 2 VT^t 



= G 



N ■ 



(16) 



i.e. we reproduce the well known result that the zero-bias conductance of the diffusive contact Gns = Gn, in contrast 
to the ballistic case when Gns = 2Gn, first obtained by Artemenko,Volkov and Zaitsevtj by a different method. 

Fig. 4 shows the results of numerical calculations of the dependence of the zero-bias conductance of a disordered 
FS contact vs spin polarization. It is seen by comparison of Figs. 3 and 4, that assuming ballistic transport in a FS 
contact one can overestimate the spin polarization in a ferromagnet. The results presented here correspond however 
to the strong scattering regime. For a more quantitative comparison with experiments the model should be further 
extended to the regime of arbitrary scattering strength. 



IV. EXCESS RESISTANCE 



So far we have taken into account both the interface and impurity scattering in the contact, but neglected the 
contribution of an electric field penetrating a superconductor. The latter can be indeed neglected in a point contact 
geometry, while it becomes important in planar contacts, in particular close to T c , when an_electric field penetrates 
into a superconductor over the macroscopically large charge-imbalance relaxation length AqEj. 

The correspondiag-eontribution to the boundary resistance of a FS contact can be calculated by the generalization 
of the approach oltl3 : E3 valid for a clean superconductor. Excess resistance R ex is given 



Rex = FXqp s /S. 



(17) 



Here p s is the normal state resistivity of a superconductor and F = Y* /Y, where Y* represents the charge current in 
FS contact 



x , f d 2 k« f +0 ° 



and Y represents the total current 



Y 



fAl [ +0 ° ( df 



A U (E)-B U (E) 



dE. 



(18) 



(19) 



Here N B (E) = Ej \JE 2 — A 2 (T) is the density of states in a superconductor. 
At E < A(T) the coefficients C, D vanish, while A, B are given by 



2M 



A hl (E) = 1 - B hl (E) = - 



T t T x (1 + q- 2 ) 



( 1 - r T r l) 2 + a2 ( 1 + r T r l) S 



(20) 
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where a = ^A 2 (T) -E 2 /E. At E > A(T) 



fc n,T 



[l + /3-r Tn (l-/?)] 2 ' 



(21) 



[l + /3-r T r;(l- 



0)] 



2 ' 



2(l + /3)(T T +r i ) 



k F1,l 



A U (E) + B U (E) + C U {E) + D U (E) = 1 



(22) 



[l + /3-r m (l-/3)] 



2 ' 



where f3 = ^ E 2 - A 2 (T)/E. 



The results of the calculations of the excess resistance factor F = Y* /Y for a FS contact are shown in Fig. 5. It is 
seen that F increases strongly at temperatui^s-jdase to T c . Given the fact that the charge-imbalance relaxation length 
Aq becomes macroscopically large near T c Ejo'EiJ we conclude that measuring the excess resistance in a FS contact 
can provide a sensitive probe for measuring spin polarization. 

In conclusion, we have presented the results of a theoretical study of interface resistance in ferromag- 
net /superconductor junctions. The Andreev reflection theory is extended in order to take into account the impurity 
scattering within the contact in the regime of strong disorder. The model is applied to the calculation of the excess re- 
sistance of a FS contact caused by penetration of an electric field into a superconductor. The latter contribution could 
be important in contacts with planar geometry and provides an additional method for measuring spin polarization in 
ferromagnets. 
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FIG. 1. Low temperature spin and charge currents in a ballistic FS contact for different spin polarizations in a ferromagnet 
for the barrier strength parameter Z=0. 
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FIG. 2. Low temperature spin and charge currents in a ballistic FS contact for different spin polarizations in a ferromagnet 
for the barrier strength parameter Z=l. 
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FIG. 3. Zero-bias conductance of a ballistic FS contact as a function of the spin polarization at various barrier strengths. 



9 



FIG. 4. Zero-bias conductance of a disordered FS contact as a function of the spin polarization at various barrier strength. 
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FIG. 5. Temperature dependence of the excess resistance in a FS contact. 
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